Deconfinement in a 2D optical lattice of coupled ID boson systems 
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We show that a two-dimensional (2D) array of ID interacting boson tubes has a deconfinement 
transition between a ID Mott insulator and a 3D superfluid for commensurate fillings and a dimen- 
sional crossover for the incommensurate case. We determine the phase diagram and excitations of 
this system and discuss the consequences for Bose condensates loaded in 2D optical lattices. 
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Loading optical lattices with ultracold atoms is pro- 
viding new types of many-body systems as well as new 
ways to look at hard problems such as quantum phase 
transitions in strongly interacting systems. The recent 
observation of the superfluid-Mott insulator transition Q 
has demonstrated the large degree of tunability offered 
by these setups. This will allow the exploration of e.g. 
spin-1 atoms, fermions, and Bose-Fermi mixtures, where 
strong interactions may lead to novel quantum many- 
body states. They offer also the possibility to realize 
anisotropic traps and thus to explore the physics of in- 
teracting particles in reduced dimensions 0J3 • Recently, 
two dimensional (2D) optical lattices have been realized, 
where at each lattice site, a ID tube of ultracold atoms is 
trapped 0, Q . Such systems thus provide a unique way 
to obtain strongly interacting ID Bose gases, of which 
the Tonks gas is the best known example (see e.g. 

In ID, interactions cause more drastic effects than 
in higher dimensions and lead, both for fermions and 
bosons, to a very peculiar state known as Tomonaga- 
Luttingcr liquid (TLL) [3, la, li| , which by now is theoret- 
ically well understood. However, despite the intense the- 
oretical and experimental activity in fermionic systems 
or spin chains |9(, much less is known when a collection 
of such ID systems is coupled. In general, one can ex- 
pect a dimensional crossover where the system goes from 
a TLL at high temperature to an anisotropic 3D system 
with more conventional properties (e.g. a Fermi liquid for 
fermions, a superfluid for bosons) at low temperature. 

The transition is even more drastic and less well under- 
stood when, due to the existence of an axial periodic po- 
tential (thereafter called "Mott potential" ) , the ID sys- 
tem becomes a Mott insulator (MI) [13, Ell E3 

Quite 

generally, the competition between Mott localization and 
superfluidity leads to interesting phenomena 01 . For 
coupled ID chains the competition between the Mott 
insulating physics and interchain tunneling changes the 
dimensional crossover into a deconfinement transition 
where the system goes from a ID insulator towards a 
3D metal or superfluid. For fermions this transition re- 
mains a theoretical challenge 0] and has been observed 
in coupled chain systems such as the organic conductors 



For bosons, much less is known even though the 
problem has been tackled in the context of spin chains 
(i.e. hard core bosons) U| or for a bosonic ladder [l7| . 
In general, however, comparison between theory and ex- 
periments is rendered difficult by the complex nature of 
the solid-state materials. 

In this paper, motivated by the clean and tunable re- 
alization of coupled ID gases offered by 2D optical lat- 
tices 0,O|, we investigate their phases and phase diagram 
(see Fig. 0. For the lattice without the axial Mott po- 
tential, we show that the system exhibits, as a function 
of temperature, a dimensional crossover from a strongly 
anisotropic normal gas to a 3D Bose-Einstein condensate 
(BEC) or superfluid, which can have a large quantum 
depletion. For this 3D superfluid (SF), we obtain the 
condensation temperature, the zero-temperature conden- 
sate fraction, the excitation spectra, and the momentum 
distribution in the thermodynamic limit. When an axial 
Mott potential is present, we show that its tendency to 
lead to a Mott insulator competes with the Josephson 
coupling between the ID tubes, which tries to delocalize 
the atoms. Changing the hopping amplitude between the 
tubes drives the deconfinement transition, and we deter- 
mine its properties and thephase diagram. In addition, 
since in the experiments [j, |5j the 2D optical lattices 
are confined in a harmonic trap, we have considered the 
effect of finite-size tubes. For small enough hopping a 
2D lattice behaves as a 2D array of Josephson junctions 
which can undergo a transition to a 2D Mott insulating 
state where tunnelling between tubes is suppressed by 
the "charging energy" of each tube. 

We consider a 2D square lattice of ID tubes of length 
L containing Nq = pqL bosons with Hamiltonian: 
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v E'r(^) is the bosonic field operator at x (the axial co- 
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FIG. 1: Zero-temperature phase diagram of a 2D lattice of 
coupled ID boson systems. The values on the vertical axis 
are defined up to a factor of order unity. K is the TLL pa- 
rameter. The corresponding values of 7 = Mg/h 2 po are also 
given. A periodic potential of amplitude uo, commensurate 
with the density po, is applied longitudinally (mo ~ 0.03 p for 
K = 1 and Uo ~ 0.02 p for K = 2.5). For infinite ID systems 
only two phases exist: a ID Mott insulator (ID MI) and a 
3D (but anisotropic) superfluid (SF). For finite tubes, a 2D 
Mott insulator (2D MI) can exist for J/p below the horizontal 
dashed curve (we have taken No = 100 atoms per ID tube). 
The vertical dashed line (schematic) indicates that the tran- 
sition to the ID MI should take place for K ~ 2. The insert is 
the phase diagram of a 2D optical lattice of finite ID systems 
for wo = and iVo = 100, as a function of chemical potential 
p and the hopping J, and for various values of K: K = 1 
(continuous curve), K = 2 (dotted curve), K = 4 (7 = 0.71, 
dashed-dotted curve). 



ordinate), and at the lattice site R = (m,ri)b, where 
m,n are integers and b the lattice parameter; pn(x) = 
^ r (x)^r(x) the density operator and u(x) is the Mott 
potential. For cold bosonic atoms confined to ID, the 
interaction v(x) — q- m Si x), where gm = 2fi, 2 a3£)[(l — 
Ca 3D /V^^±)M£ 2 ± ]- 1 [3, being a 3D the 3D scattering 

length, C ~ 1.4603, and t±_ ~ bV Q ^ /4 /ir the trans- 
verse oscillator length Q (Vox is the strength of the 
transverse optical potential in units of the recoil energy 
E R = h 2 it 2 /2Mb 2 ). In the last term of 0, which de- 
scribes the Josephson coupling of the ID tubes, (R, R') 
stands for sum over nearest neighbors. Furthermore, we 
assume that the 2D optical lattice is deep enough for 
the hopping J -C p,, p, being the ID chemical poten- 
tial. We first focus on the thermodynamic properties 
(i.e. No — > 00, and L — > 00); finite-size and trap effects 
will be considered towards the end. 

For the low-temperature properties of lJT]l. it is con- 
venient to use the so-called bosonization technique 0, 



I @. Using * R (a) ~ V^Pu" e^ R(a:) and p R (x) = 
[po + d x 8 R (x)/ir] J2 m e 2lm ^ R ^ +7rpoX l Eq. Q becomes 
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where <5 = 2po — G/n is the mismatch between the den- 
sity po and the periodicity of u(x) = uoCOs(Gx). The 
dimensionless couplings gj — 2irJ(poa) 2 A/hv s po, and 
g u = 2iruo(poa<) 2 /hv s po, where a « hvJfi is the short- 
distance cut-off and A « (Kjn) 1 ! 211 1\. Results and 
asymptotic expressions for K and the sound velocity v s 
in terms of 7 = Mguj/h 2 po can be found in Ref. 0. In- 
terestingly, an anisotropic version of the Bose-Hubbard 
model Il3j (i.e. with axial hopping J x 3> J) also leads to 
Eq. J2J) |l9j , but the relationship of the microscopic pa- 
rameters (namely J x , J and U) to K, v s , g u and gj is not 
easy to estimate. But if instead K and g u are regarded as 
phenomenological parameters, our results directly apply 
to this model as well. 

First we consider the case uq = (i.e. g u = 0). This 
includes the incommensurate case <5 ^ for which the 
periodic potential u(x) has no effect at low tempera- 
tures 1 1 CX 1 1 lL 1 1 2| . In the absence of tunnelling, the system 
is then a set of isolated TLL's with no true condensate 
and a power-law decay of phase correlations |lfj| |. Tun- 
neling thus induces, at a temperature T c , a dimensional 
crossover from the TLL behavior at T > T c , towards a 
3D superfluid or BEC at low temperatures. To study 
the crossover, we treat the tunneling term in a mean- 
field (MF) approximation, taking i/)o = (^r(#)) to be 
the (square root of the) condensate fraction. The prob- 
lcm then reduces to a sine- Gordon (SG) model [lflllii, 
describing an effective ID system. Above T c , ipo = 0. 
Therefore, T c is found by the existence of a non-zero so- 
lution for ipo. We get, in agreement with [2f1 |. 
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where f(K) = Asm (^) B 2 (^, 1 - and B(x) is 
the Beta function. Below T c , ipo is found by minimizing 
the free energy energy density. At T — one minimizes 
£ M F(^o) = 4J^-^^tan. 



where the second 



ihv s 2(&K- 

term is the SG model energy density [21|, and A s (t/j ) the 
soliton gap (the excitations of the SG model are gapped 
solitons and breathers H H 0). Thus we find that 
the condensate fraction ip 2 (T = 0) ~ p ( J/p) 1 /^-^ 
whereas the soliton gap A s - p(J / p) 2K / {AK |2sj . 

In the broken-symmetry (superfluid) phase the system 
possesses two modes. The Goldstone mode (i.e. oscilla- 
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tions of the order parameter phase) is gapless as a con- 
sequence of (global) gauge invariance. The other mode 
corresponds to oscillations of the order parameter ampli- 
tude. To obtain these modes we compute the gaussian 
fluctuations around the MF solution (random-phase ap- 
proximation, RPA). Using known properties of correla- 
tion functions of the SG model 0,0 Eli an d within the 
single-mode approximation (SMA: taking into account 
only the lowest breather mode) 0, we obtain: 



w+(?,Q) 



A 2 , 2 2, 

A+ +v s q + 



A 2 + /2h 2 
8K-2 



F(Q) 



(4) 
(5) 



where A 2 + = A 2 (8K - 2)/(8K - 1) and q (rcsp. Q) 
is the momentum along (resp. perpendicular to) the 
chains. A„ = A s sin (mr/2(8K — 1)) is the energy gap 
for the n-th breather (n — 1,2,...), and F(Q) — 
^ (1 — cosQjb). 
j=y,z 

For weakly interacting bosons (K ^> 1), the use of 
SMA may be questioned because the breather excita- 
tions of the mean-field SG model proliferate @, 0,113 • 
Therefore, we have employed a variational approach pre- 
viously used for coupled spin chains 23], and found the 
same dependence on J/p, for the dispersion as in 
With this approach, we have also calculated the momen- 
tum distribution of the superfluid phase (in the absence 
of a trap). At T — and small momenta (<gC 7r/6), 
* R (x) ~ pl /2 e^ R W and u£(c?,Q) ~ v 2 q 2 + v 2 ± Q 2 , 

(v ± ~ ixKJIv) 2K,{iK ~ l) /K, cf. Eq. Q), we find an 
anisotropic generalization of the Bogoliubov result : 
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which is valid for arbitrary interactions within the tubes. 
w(Q) is the Fourier transform of the Wannier orbital, 
which varies slowly with Q. 

We next consider the case of a finite axial Mott poten- 
tial, i.e. u 7^ and S ~ 0. In the absence of tunneling, 
the individual tubes would be ID Mi's: even for a weak 
Mott potential, an integer number of bosons are localized 
to each well of the Mott potential 0, 0, ^3 • Increas- 
ing the tunneling will thus lead to deconfinement (i.e. 
ID MI to 3D SF transition). To study the competition 
between the the Mott potential (g u ) and the Josephson 
coupling (<?j), we have used the renormalization group 
(RG) method. Since the Hamiltonian in defines an 
effective field theory, its couplings (K, g u , gj, etc.) de- 
pend on the cut-off energy scale, which is typically set by 
the temperature for T < /j,. Thus, as the temperature is 
decreased, the couplings K,g u ,gj are renormalized, due 
to virtual transitions to high energy states. The nature 



of the ground state is thus determined by the dominant 
coupling as T — > 0. Using perturbative RG, we find the 
following RG flow equations: 
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where I w ln^/T. The coupling gp is generated by the 
RG and describes an interaction between bosons in neigh- 
boring tubes; g u and gj compete for 1/4 < K < 2: the 
relative magnitude of their initial (i.e. bare) values de- 
termines which coupling grows faster. The faster growth 
of one coupling inhibits the other's growth via the renor- 
malization of K: If g u grows faster, K flows towards 0, 
leading to the ID MI, while if gj grows faster, K also 
grows and leads to the 3D SF. To estimate the phase 
boundary between the ID MI and the 3D SF, we inte- 
grate these equations (at fixed g u {0)) with various g.j(0) 
till the couplings are of order one. The resulting phase 
boundary is the continuous curve in Fig. ^ Via a sep- 
arate mean field calculation at a specific K value (see 
(l^l. we showed that the condensate fraction Vo grows 
continuously from zero at the phase boundary. 

Finally we turn to finite-size effects. First the finite ex- 
tent of trapped cloud (either longitudinally or transver- 
sally) limits the minimum momenta of the modes. Thus 
the energy of the lowest modes in the 3D SF phase can 
be directly estimated from our results by using the min- 
imum available momentum in (I4I5II . For instance, for 
a finite 2D lattice containing M y x M z tubes (i.e. an 
atom cloud of size L x M y b x M z b) the lowest available 
momentum is ~ w/(Mib). Putting this value into (1 ll- r >[) 
shows that the frequency of the lowest transverse modes 
decreases with decreasing J. This is in qualitative agree- 
ment with the hydrodynamic analysis of Ref. |25| . How- 
ever, unlike this hydrodynamic analysis, our results apply 
to the strongly interacting regime of the superfluid near 
the transition where depletion of the condensate can be 
large (as shown e.g. by the power-law behavior of the 
condensate fraction, i\) 2 ~ po(J/ /J,) 1 '^^ 1 '). In addi- 
tion, we are also able to describe the transition to the 
Mott insulating regime. Note that for finite size systems 
the phase transitions described above become crossovers. 
Furthermore, for harmonically trapped systems in the MI 
regimes an inhomogeneous state arises where the the SF 
and MI coexist 0,E3i an d V'o decreases slowly across the 
phase boundaries. 

A second type of finite size effects comes from the finite 
size and hence discrete spectrum of each tube. In the 
experiments in 2D optical lattices loaded with ultracold 
atoms 0, 0, there are typically N ~ 10 2 — 10 3 atoms 
per tube. At low temperatures, if J is made very small 
(i.e. the lattice potential very deep), each ID tube will 
behave as an atomic "quantum dot" characterized by a 
charging energy Ec = hirv s /KL 0. This is to be 
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balanced against the hopping energy Ej, which is not 
just J No as the hopping amplitude is renormalized by 
the creation of virtual sound waves within the ID tubes. 
Thus, to study this effect for u = 0, we integrate out the 
phonons perturbatively in J and obtain a quantum-phase 
model: 

H QP = -Ej cos (^0R - </>0R') 

<R,R'> 

R R 

where ATr (resp. 4>or) is the particle-number (reps, 
phase) operator of tube R (hence [Ar,0or'] = «<5r.R')- 
The renormalized hopping Ej = Ej(Nq) ~ JA 2J U 
The model © has been extensively studied in the litera- 
ture in connection with 2D arrays of Josephson junctions. 
It exhibits a 2D MI-SF at a critical value of Ej/Ec- For 
commensurate filling No (fi — 0), a Monte Carlo (MC) 
calculation j2|| gives (Ej/Ec) c — 0.15. In our case this 
reduces to ( J//i) c — 0.3A 3 ^ 2 in the Tonks limit (K ~ 1) 
and to (J//i) c — 0.15Aq~ 2 for weakly interacting bosons 
(K ^S> 1). Away from commensurate filling (/i ^ 0), the 
critical J//x is reduced as shown in the insert of Fig. ^ 
obtained using the MC data from Ref. |26j. Note that 
in the 2D MI, only the phase coherence between different 
tubes of the 2D lattice is lost. There is no gap (apart from 
the finite size gap) to excitations within each ID system. 
This is different from the ID MI described above, where 
all excitations are gapped and any order (even quasi- 
long range) in the phase is absent. For finite uq, when 
K ^ 2, g u flows to zero and the consideration above ap- 
plies (hence the horizontal and vertical dashed lines in 
Fig-Q. Otherwise [K < 2), gj flows to zero and the ID 
MI is as described above. 

The above predictions can be directly tested in 2D op- 
tical lattices |4j, |5j. In presence of the Mott potential, 
the deconfinement transition is of course the easiest to 
check. Even without the Mott potential, V'o an d T c as 
functions of J/V, and the change in the momentum dis- 
tribution from the ID TLL to the 3D superfluid should 
be accessible, by releasing the trap and the optical lat- 
tice, then measuring the expansion images. Finally, the 
excitation spectra Q4I5JI of the superfluid can be measured 
with Bragg spectroscopy at low momentum transfer. 
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Note Added.- After this Letter was submitted, we re- 
ceived a preprint from Esslinger reporting on exper- 
imental results for the system we have considered here. 



The reported results are at least qualitatively consistent 
with our predictions. 
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